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EDGE PERCOLATION ON A RANDOM REGULAR 
GRAPH OF LOW DEGREE^ 

By Boris Pittel 
Ohio State University 

Dedicated to Jimmy and Katie Lyons 

Consider a uniformly random regular graph of a fixed degree 
d > 3, with n vertices. Suppose that each edge is open (closed), with 
probability p{q — 1 — p), respectively. In 2004 Alon, Benjamini and 
Stacey proved that p* = {d — is the threshold probability for 
emergence of a giant component in the subgraph formed by the open 
edges. In this paper we show that the transition window around 
p* has width roughly of order n~^^^. More precisely, suppose that 
p = p{n) is such that cj := n^^^ \p — p*\ — > oo. If p < p* , then with high 
probability (whp) the largest component has 0((p — log n) ver- 

tices. If p>p* , and logLj 3> log log 71, then whp the largest compo- 
nent has about n(l — {pn + g)'*) x n{p — p*) vertices, and the second 
largest component is of size (p — p*)^^(logn)^^°'^\ at most, where 
TT — {pn + q)''^^,TT G (0, 1). If Lo is merely polylogarithmic in n, then 
whp the largest component contains 1^^^^^+°^ vertices. 

1. Introduction and results. Let d>3 and n > d be given. Assuming 
that nd is even, introduce Gn, tiie sample space of all d-regular graphs on 
the vertex set [n]. It is known (Bender and Canfield [5]) that 

,^ , M-1)!! f 1, 
(1.1) \Gn\ ='-J^^^P[ — + Oin~') 

Here is a probabilistic interpretation of (1.1), due to Bollobas [7], who ac- 
tually considered a general degree sequence. Introduce n sets Si,. . . ,Sn, Si 
representing vertex i, with \Si \ = d. Then (nd — 1)!! is the total number of 
complete matchings (pairings) on S = [j^Si, and the exponential factor is 
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the probability that the matching chosen uniformly at random is graph- 
induced, that is, there are no pairs (s,t) with s and t from the same set Si, 
i G [n], and no pairs with s,s' G Si, t,t' £ Sj, i / j. Equation 

(1.1) follows then from the observation that every d-regular simple graph 
induces (d!)" distinct pairings on S. 

Consider Gn the uniformly random ti-regular graph on n vertices. Thus 
Gn assumes every one of its \Qn\ values with the same probability |^n|~^- 

Let p G (0, 1) be fixed. We define the random graph Gnp as the random 
subgraph of Gn obtained by "opening" ("closing") each edge of G„ with 
probability p ((? = 1 — p), independently of all other edges. Thus the edge set 
of Gnp is the set of all open edges of G„. Several years ago Itai Benjamini 
[6] posed a problem (i) to show that p* = {d— 1)"^ is a threshold value of p 
for emergence of a giant component in Gnp, and (ii) to determine the width 
of the transition window around p* . In 2004 the part (i) was solved by Alon, 
Benjamini and Stacey [2]. 

Here is a quick-and-dirty argument for why p* had better be the threshold. 
A subgraph of G„ induced by the vertices within a relatively small distance 
from any given vertex u is a tree in which all nonleaves have degree d. 
Orienting the edges away from v, we get a directed tree in which the root 
V has outdegree d, and the remaining nonleaves each have outdegree d — 1. 
And it is known (Durrett [13]; see also Grimmett [17]) that, for the edge 
(bond) percolation on an infinite directed tree of outdegree d — 1, the critical 
probability is p* = {d— 1)^^, and for p> p* the probability that the root is 
in an infinite cluster is 1 — vr, 



Here vr = 7r(p) is the probability of eventual extinction for the branching 
process with immediate family size having binomial distribution Bin(d — 
l,p). Thus we should anticipate that each of the d neighbors of vertex v 
will have only few "descendants" with probability tt, independently of other 
neighbors, and so v itself will have only few descendants with probability 



Guessing that with high probability (whp), that is, with probability 1 — 
o(l), each vertex v is either in a small component or in a unique "giant" 
component, we can even anticipate then that the fraction of vertices in the 
giant component is asymptotic to 



(1.2) 



vr = {np + q) 



d-l 



^G(0,1). 




(1.3) 



a{p) = 1 - (pvr + q)'^. 



(This was not proved in [2].) 
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What about the window width? Since the pioneering work of Erdos and 
Renyi [16] in the early 1960s, it has been known that for d = n — 1, that 
is, Gn = Kn, a unique giant component appears when the average vertex 
degree d exceeds 1. The issue of the transition window around 1 remained 
wide open until Bollobas [9] (see also [10]) was able to show that its width 
is of order 0(n~^/^(logn)^/^). Later Luczak [19] extended the techniques of 
[9, 10] and showed, among other sharp results, that the logarithmic factor 
could be removed. (See also [18, 21, 25] concerning the critical phase d — l = 
0{n~^^^).) In light of these results Itai Benjamini wondered whether the 
bound n~^/^ would hold for the graph Gn of a small fixed vertex degree 
d as well. Recently Borgs et al. [11, 12] have established this width bound 
for a class of deterministic regular graphs with vertex-transitivity property, 
meeting a certain "triangle condition." This class contains Kn and some 
"high-dimensional" tori. Notice that whp the random d-regular graph G„ is 
not in this class, since almost all d-regular graphs are asymmetric (Bollobas 
[8]), whence intransitive. 

Our goal is to solve part (ii) of Benjamini's problem, confirming his n~^^^ 
conjecture. Here are our main results. 

Given p, let Ln\Ln^ denote the size of the largest component and the 
size of the second largest component in Gnp- 

Theorem 1 (Subcritical case). Suppose thatp = p{n) is such that 

(1.5) lim n"'^/'^(n* — p) = oo. 

n^oo 

Then (^p*_py^2i^gn bounded in probability, or Ln^ = Op{{p* — p)~^log?i) 
in short. 

Theorem 2 (Nearcritical case). Suppose that, for some a> 0, 

(1.6) n^/^\p-p*\<(logn)''. 
Then, for any b > max{a, 1/3}, 

(1.7) lim P(?i2/3(iogn)-2'^ < < n^/^(logn)'') = 1. 

n — ^oo 

Theorem 3 (Supercritical case). Suppose that 

(1.8) ni/=^(p-p*) > (logn)"("), 
where a{n) oo however slowly. Then 

= (1 + op(l))n(l - [pTT + q)'') n{p - p*), 

(1.9) 

4'^ = (p-P*)-'(logn)i+°^W; 
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here op(l) denotes a random variable converging to in probability, and the 

second equality for Ln^ means that, in probability, Ln^ is of order n{p — p*) 
exactly. 

Asaf Nachmias and Yuval Peres (personal communication) informed me 
that they have obtained similar (but somewhat more precise) results. 

These estimates are very similar to, but not as sharp as, those in Bollobas 
[9, 10] and Luczak [19] for G„ = Kn- In both papers the proofs relied heavily 
on Bollobas' striking bound 

(1.10) C{u, iy + l)< (c/£)^/2z.'^+(3^-i)/2, 

C{v,v -\- 1) being the total number of connected graphs with u vertices and 
u + l edges. (See Luczak [20] for a short proof of (1.10) with the best constant 
c.) While there is a formula for the number of trees with a given degree se- 
quence (Moon [24]) no analogue of the general bound (1-10) for the number 
of connected graphs with a given degree sequence seems to be known. Fortu- 
nately it is possible to obtain sufficiently sharp estimates for the parameters 
of a randomized algorithm that determines a graph component containing a 
given vertex. These estimates, together with an asymptotic formula for the 
likely count of certain tree components in the case p = p* , yield the proba- 
bilistic bounds for Ln^ and Ln \ More specifically, we derive the equations 
for the conditional (one-step) expectations of this process, which suggest 
a system of deterministic differential equations whose solution should be a 
good approximation for the scaled random trajectory. That this is indeed 
the case is proved by using exponential (super)martingales constructed from 
the integrals of these differential equations. More general, martingale-based, 
techniques for random graphs developed by Wormald [27] do not provide 
estimates strong enough for our purposes. Our method is reminiscent of 
a technique used in [26] (A;-core problem), [3] (Karp-Sipser matching algo- 
rithm), and [1] (random graphs with immigrating vertices). But most closely 
related is our joint paper with Jozsi Balogh [4] in which we used the dif- 
ferential equations and the exponential supermartingales for analysis of a 
"bootstrap" site percolation on the random d-regular graph. We must also 
mention that Molloy and Reed [22, 23] had used the differential equations 
in their study of existence of a giant component in a uniformly random 
graph with a given degree sequence. A more or less direct application of 
their remarkable result to Gnp seems fraught with subtle difficulties, even 
for p bounded away from p* . Certainly our bounds for the transition window 
width could not be obtained this way. However, an algorithmic approach to 
the giant component phenomenon ushered in [22, 23] is a key tool in both 
[4] and the present paper. 

We believe that our techniques can be extended to a more general (uni- 
formly) random graph Gn,d with a given degree sequence d = (di, . . . , dn), 
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such that each di G [3,d], d < oo. Intuitively, in this case the counterpart of 
(1.2) is 

j=3 l^i'e[n] "i' 

(can the reader see why?), and so 

1 

^ E,jA,-i- 

And (1.3) becomes 

j 

We think that again the transition window around p* has width of order 
n-V3. 

In conclusion we should mention a new book by Durrett [15] in which the 
random graph processes are used systematically either as proof tools or as 
probabilistic models of evolving networks. The percolation on Gn^d we have 
just described is inspired by the discussion of the largest cluster problem for 
a random graph, with a given degree distribution, in Chapter 1 of this book. 

The rest of the paper is organized as follows. In Section 2 (p close to p*) 
we evaluate sharply two first moments of the number of tree components of 
sizes dependent on \p — p*\, and deduce the lower bound for Ln \ that is, 
the left half of (1.7) in Theorem 2. In Section 3 we introduce the percolation 
process on the random pairing and show that it is enough to analyze this 
eminently more tractable process. Motivated by the formulas for the condi- 
tional expected state changes of the attendant Markov chain, we introduce 
and solve the deterministic system of differential equations, and prove — via 
a family of the exponential supermartingales — that whp the random process 
stays sufficiently close to the deterministic trajectory. In Sections 5, 6 and 7 
we use the probability tail estimates to complete the proof of the Theorem 2, 
and to prove Theorems 1 and 3. 

2. Proof of Theorem 2 (Nearcritical case). Lower bound. We begin at 
the middle since the proof is purely enumerational, a natural extension of 
an argument for Kn by Bollobas [9, 10]. 

First of all, from the discussion of (1.1) in the Introduction we know that, 
for the uniformly random pairing Vn on the set S = IJiLi Si, and the event 
An = {Vn is simple graph- induced}, 

(2.1) P(^„)=exp(^-^^ + 0(n-i)^ ^ exp (^-^^) > 0. 
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So we will study percolation on the random pairing Vn instead of percola- 
tion on the graph Gn- Here we open each pair {s,t) G Vn with probability 
p, independently of other pairs. And the problem becomes: find the proba- 
bilistic bounds for the sizes of the largest component and the second largest 
component of the (multi)graph MGnp on [n] determined by the subpairing 
Vnp consisting of the open pairs. [Each {s,t) G Vn, s G Si, t £ Sj, gives rise 
to an edge if j, a loop at z if i = j-] According to (2.1), for every 
set 7i of graphs on [n] , 

(2.2) F{Gnp G W) = FiMGnp G n\An) = 0(P(MG„p G Ti)). 

Thus any rare event for the random multigraph MGnp is a rare event for 
the random graph Gnp- 
Let 

n^/^\p-p*\ < (log?i)". 
We need to show that whp 

41)>7l2/3(logn)-2-. 

It is enough to prove that whp the size of the largest tree component exceeds 
77,2/3 (log n)"^"^. 

For A: > 1, let E'fc denote the expected number of tree components of MGnp 
with k vertices, and for k^,k'^ >1 let Ef.i ^,2 denote the expected number of 
ordered pairs of tree components, with and k'^ vertices, respectively. 

Lemma 1. (i) If k = o{n'^/^), then 

(2.3) . „./-.,-.-..^^^M^e-. 

(2.4) ^„^:=a— +o(l), 

n 

(2.5) a:= — (p-p), P = j-^- 

(ii) If k^,k^ = o{n^/''^), then 

(2.6) |^ = l + af^-E^)+o(l), k = k^ + k^. 

Notes. (1) Observe that in (2.5) a = for p = p*. (2) For k = 1, we 
obtain 

So the number of isolated vertices, with no loops, is whp asymptotic to nq'^, 
which should be expected, of course. 
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Proof of Lemma 1. (i) By symmetry, 



(2.7) Ek=[ljPk, 

where is the probabihty that [A;] = {1, . . . , A;} is a vertex set of a tree 
component. Let us compute this probabihty. Let S = {6i, . . . ,6k) denote a 
generic degree sequence of a tree on [k]. It is weh known that 5 > is a 
degree sequence of such a tree iff 

je[fc] 

and that the number of trees with a given feasible S is [24] 

(fc-2)! 

Besides, each 6i is at most d, since we are deahng with a subgraph of Gn- 
For a corresponding subpairing of Vn we need to pick 6i points from the set 
Si, i £ [k]. Once it is done, the number of tree- induced subpairings on the 
chosen points is 

Y^^^^^^ ■ n S.l = ik- 2)1 Yl 6.. 

n.e[fc](5.-l)! ^m 

Let Ai,...,Afc be the generic numbers of additional points from the sets 
Si,...,Sk which are paired among themselves rather than with points in 
the sets Sk+i, ■ ■ ■ , It is necessary, of course, that 5j + Aj < d, i e [k], and 
that A := J2i£[k] '^i even. The number of ways to pair these A points is 
(A — 1)!!. The remaining di = d — 6i — Ai points in each set Si, d = J2ie[k] 
in total, are to be paired with some d points from Ui^[fc] ^i, while the rest 
of points from Ut.^[fe] Si, all (n — k)d — 9 of them, are to be paired among 
themselves. Combined, this can be done in 



nd — kd 
d 



^ d\{{n - k)d - d - 1)\\ 



ways. In addition, the number of ways to split each Si into three ordered 
subsets of cardinalities 6i,Ai,di is 



Si,Ai,di 



The probability that a tree subgraph, with the degree parameters S,A.,d, 
is a component of the random sub (multi) graph MGnp is 
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Let P{S, A.,d) denote the probability that Gnp has a tree component with 
vertex set [A;], and the degree parameters S,A.,d. Collecting all the pieces, 
we obtain 

^^■^^ = (nd-iyjS '' ~ ^^-^^ " ( ^•^^'^ -k)d-d- 1)!! 

X n 6id5i,Ai,di-p''-'q^/^+^. 

iG[fe] 

Notice that A uniquely determines d, 

d = kd-{k-l)-A. 

Let -P(A) stand for the probability of a tree component with a given A; 
so -P(A) is the sum of P((5,A,9)'s over (5's, A's and 9's satisfying the 
conditions 

Ai + Si + di = d, ie[k], 
(2.9) ^5, = 2(A;-1), 

^A, = A. 

ie[k] 

According to (2.8), given A, we need to evaluate 

meet (2.9) 

= [-^'1 E n 

aAa>0 ie[fc] ^*■ 

ai+^5i+9^=c^-l,^G[A,•] 

ai+Ai+ai=d-l,iG[fc] 

^ [xV]((l + x + y)('^-i))'- 



1 (A:(d-1))! 



[(d-l)!]'= (fc-2)!A!5!' 
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Therefore, summing P(d,A,9)'s over all S's, A's, d s for a given A, we 
obtain 

P,^) . ^:W^-'>>"^:''" ( -/'^l ((n - - a - 1)11 



(nd-l)!! A! V 5 

(2.10) X 

9 = A;d-2(A;- 1) - A. 

Using (2s — l)!!/(2s)! = (2^*5!)"^, substituting A = 2D and excluding d, we 
rewrite (2.10) as follows: 

^^^^^_ ^fc2fe'^-^'+Hfc(d-l))!((n-fc)d)! .,M-2k+2 
2"rf/2(nd-l)!! 

(2.11) 

(4<?)-^ 



D\{kd -2k + 2- 2D)\{nd/2 -kd + k-l + D)]' 

What is left is to find a sharp estimate for the sum of P{2D) over D < 
kd — 2{k — 1). Notice that the D-dependent fraction in (2.11), call it F{D), 
is bounded above by 

^+.^^ I 1/ (fcd-2(fc-l))^ \^ 

^ ^ {kd - 2k + 2)l{nd/2-kd + k- 1)1 D\\Aq{nd/2 - kd + k - 1) J ' 
Introduce a Poisson(Afc) random variable Z, where 

, , (kd-2(k-l))^ k'^{d-2f , 

Observe that is of order A:^/n exactly. Recall that k^ jr? 0. Using 
Chebyshev's inequality, 

,2\ e\ / 2e 



P(z....-(S)>(5) 



n ^ oo, 



for every fixed e > 0. Therefore, denoting Afc + X]/'^ {n'^ /k^Y by Dfc, 
(2.13) E^^(^)< 



Z)>D, 



fe 



(A;(i-2/c + 2)!(n(i/2-A;d + /c-l)! Vn 



.2 



For Z) < Dfc, we have 



2D 

[{kd -2k + 2- 2D)l]^^ = [{kd -2k + 2k)\]-^ Y[{kd-2k + 2- j) 

12 , 



[{kd -2k + 2)\]-\kd -2k + 2f^e°^^ll^\ 
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where 



0, 



k k k \k'^ J \n'^ xn"^/^ , 

provided that e = 1/6, say. So 

[{kd -2k + 2- 2Dy.]-^ = (1 + o{l))[{kd -2k + 2)!]~^(A;d -2k + 2f^, 
and hkewise 

[{nd/2 -kd + k-l + D)l]~^ = {1 + o{l))[{nd/2 - kd + k - . 
Consequently 

y: F{D)={i+o{i)) y: F^iD) 

D<Dk D<Dk 

^'■''^ = + -W \.d-2fe + 2)Knd/2-fed + fe-l)! ^(^ < 



(1 + 0(1))- 



(/cd -2k + 2)\{nd/2 -kd + k-l)\' 
Combining (2.11), (2.13) and (2.14), we get 

D<{kd~2k+2)/2 

(2.15) 

_ d^2^'^-^+^{k{d - l))!((n - h)d)\p^-^q^'^-'^^+^ 

- ( + oU) 2nd/2(„j _ i)\\Q.d -2k + 2)\{nd/2 -kd + k-l)\'^ ' 

Call the last fraction Tj.. Tnk-, the n-dependent portion of Jtj, can be rewrit- 
ten this way: 

_ {{n-k)d)\ 

nk 



2nd/2^nd - l)!!((n/2 - k)d+{k - 1))! 

{n-k)d + k-l 

{{n-k)d)\ \ ^ 



{{n - k)d + k - 1)1 ( nd 

I nd 

\ 2 

Evaluating each of the two factors on the right in a standard way, we obtain 

^nfc = {ndy^+^2~^'^+^~^e^-\ 

where 
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Therefore the formula for Tk simpUfies to 

(2.17) = n-(>'~'^dp''~V'~"^'~'^ ^^f}'l~^ - 1))! ^^"'- 

Using (2.7), (2.15), together with (2.12), (2.16), (2.17), and 



n\ n' 
k 



we obtain then 



(2.18) E,-[j^jP,-ndp q \^^M-2{k-l))\' 



where 



^2 

*nfc = Afc + $„fc- — + 0(1) 

In 

= a— + o(l), 
n 

a := — (p -p), p 



2dq - " - 

The proof of part (i) is complete. 

Part (ii). Since the proof parallels the above argument, we will skip some 
details. First of all, 

(2.19) ii;fci,fc2= (^^i^^2'^^_^)Pfci,fc2, k:=k^+k\ 

where P^i ^2 is the probability that the sets \k}\ = {1, . . . , fc^}, [/c^] := {fc^ + 
1, . . . , fc} are the vertex sets of the tree components and in MGnp- 
Let ((5^, A^,9^) and ((5^,A^,9^) denote the degree parameters of and 
T^, as in Part (i). Then P[6, A.,d), the probability that MGnp has the tree 
components and with those parameters, is given by 

PiS,A,d) 



{nd 

5" A" d 



X 



(2.20) 



a 



a I \ a 



12 
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Vi9 -a,9^-a,(n-/c)(i-9 + 2ay^J^-^^ ^ 

X {{n-k)d-d + 2a-iy.\, 

where d = + d^. The first hne expression contains a double dose of some 
of the corresponding factors in (2.8). In particular, 

s=l 

is the probability that all the edges of and are open, and all other 
edges induced by the vertex sets [k^] and [/c^] are closed. As for the second 
and third lines sum, the ath summand is the number of ways (1) to choose 
a points among points from UiG[fci] '^i a points among points from 
UiG[fc2] Si and match them, then (2) to choose among the remaining (n — k)d 
points — a partners for the remaining 0* — a points from UiG[fc'] ^i, t = 1,2, 
and finally (3) to match the rest of the (n — k)d — d + 2a points among 
themselves, multiplied by the probability q^~'^ that all the pairs {u,v), such 
that u G UiG[A:*] Si, V ^ Uje[fe*] Si, t = \,2, are closed. 

Given numbers D*, let P{2D^ ,20"^) denote the probability that MGnp 
has these tree components and with the parameters A* = 2D'', t = 1,2. 
Summing P{d,A.,d) over feasible 5's, A's and d s, we obtain 

P{2D^,2D^) 

_ d''{{n-k)d)\2'"^-''+^Ul=i{k'{d-l))\p''-^q'"^-^''+^ 
~ 2"'^/2(n(i- 1)!! 

(2.21) xS{D^,D^), 



S{D\D^) 

r)~2D-a^,-D-a 

_ \ ^ f y 

~ a a! nLi D%ki{d -2) +2- a- 2L>*)!(nd/2 - kd + k - 2 + a + D)\' 

where D = + D^. Acting as in part (i), we upperbound the oth summand 
in (2.21) by 

1 ( Y{l=i{k'd-2k^ + 2) Y ^ 1 / {k^d-2k^ + 2f 
'^.\q{nd - 2kd + 2k - A) ) ' 'D^.\2q{nd - 2kd + 2k-4)) 

X (jj{k\d - 2) + 2)1^ {{nd/2-kd + k-2)l)-\ 
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Then the sum of the first Hne product over a, , is at most 



exp 



nLi(fc*(^-2) + 2) 



+E 



{k\d-2) + 2f 



q{nd - 2kd + 2k - 4) f-^ 2q{nd - 2kd + 2k - 4) 



- exp 



: exp 



qnd ^ 

(^(^^-2))^ 
2qnd 



2qnd 
+ 0{k/n + k^/7?) 



k = k^+k^. 



And, in fact, that sum is asymptotic to the last expression. This can be 
proven by using three independent Poissons, with parameters 

{k\d-2) + 2f 



(t = l,2), 



2q{nd-2kd+2k-A) 

ULiik\d-2) + 2) 
q{nd-2kd + 2k-A)' 

compare the proof of part (i). 
Using (2.21), we obtain then 

^ 2"<i/2(nd- l)!!nLi(fc*(d-2) + 2)! • {nd/2 - kd + k -2)\ 
■{k{d-2)r 



X exp 



2qnd 



Here 



{{n-k)d)\ 



2'^<i/^{nd-iy.l{nd/2-kd + k-2)\ 

I n 2nd 2qnd 
compare part (i). Combining the last two equations and 



+ 0{k/n + k^ /n'^^ 



n 



n 



k^,k'^,n-kj kMk^l 
we obtain 

Ef^i 1.2 = nk^ , k'^ ,n — kP^^i ^2 



exp( -— +0{k/n + k^/n'^' 



2n 



n 

t=i 



{k\d-l))\ 



kn{kH-2{k^ - 1))! 



^ak'^ /n 
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This relation and (2.3), (2.4) imply (2.6). □ 

Corollary 1. Suppose that p — p* — > 0. Let uj = uj{n) oo, such that 
u = o{n^^^), uj>n^^^\p-p*\. 

Then 

lim P(4^) > n2/3/c^^) = 1. 

For 

n^/^\p-p*\ < (logn)", 
the corollary implies that whp Ln^ is at least n^/^(logri)^^'^. 

Proof of Corollary 1. Let kn= [n^/^/w^]. Introduce X„ the to- 
tal number of tree components of MGnp of size k € [kn,2kn]- Then, by 
Lemma l(i), 

(2.22) E[X„|. E E,.„.^ 

Here, by (2.5) and the condition of Corollary 1, 

k^ 

(2.23) a— = Oiuj-^). 

n 

Furthermore, by the Stirling formula for factorials. 



{k{d-l))\ 1 



(d-l) 



d-l 



exp(0(l)). 



k\{kd-2{k-l))\ A;5/2 
if A; — > oo . Observe that 

{d-2)d-2 ^ ' J ' 

and that the function 

/(P)=M''"'=P(1-P)^"' 
attains its absolute maximum at p = p* , with f"{p*) > 0. Hence 

= exp[0(A:(p - p*f)] = exp(0(l)), 

uniformly for k € [A;„,2A;„]. Therefore E[X„], given by (2.22), is of an exact 
order 

2/Cri 
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Further, by Lemma l(ii) and (2.23) 



( 



2k, 



) 



2 



E[X„(X„-1)] 



'k\k^ 



Therefore 



var[Xj = o((E[X„])2), 



and, by Chebyshev's 



inequahty, 
hm P| 



E[Xn] 



-1 <e =1. 



Thus P(X„ > 0) ^ 1, and it remains to notice that L^!:' > n'^/^/uj'^ if X„ > 0. 
□ 

3. Percolation process on the random pairing. The upper bound for Ln^ 
in Theorem 2, and Theorems 1, 3 wiU be proved via analysis of a stochastic 
process which describes growth of a component of MGnp containing a given 
vertex, or more generally, growth of a set of vertices that can be reached 
from a given set of vertices via open edges. This growth process for MGnp 
is defined through its counterpart on the random pairing Vn', compare [4]. 
It will be convenient to assign the labels 1, . . . , nd to the points of UiG[n] '^i- 

3.1. Percolation as a Markov chain. Let ^(0) C [n] be given. We inter- 
pret ^(0) as an initial set of active vertices. Then 2^(0) = [n] \ A{0) is an 
initial set of inactive vertices. For each i € ^(0) [i £ resp.], we have the 
set 5^(0), of cardinality d, of active (inactive, resp.) points. We are about 
to define a process {A(t),I{t), Si{t),i € A{t) U I{t)}t>o, where Si{t) C Si. 
Naturally, Uje^(t) 'S'«(*)i UiGX(t) '^«(*) called the set of currently active 
points, and the set of currently inactive points, after t steps. 

At step t + 1 we (a) choose an active point s' G [jieA{t) 'S'j(Oi with the 
smallest label, say; (b) identify s" , the partner of s' in the random pairing 
Vn and delete both s' and s" from the two sets Si'{t), Si'/{t) that contain 
them; (c) if i" € I{t) and the pair (s' , s") is open, then vertex i" is pulled from 
I{t), {lit + l)=I{t) \ {i"}), and added to A{t), {A{t + 1) = A{t) U {s"}), 
so that all the points in 5j"(t + 1) := Si"{t) \ {s"} become active. 

Clearly A{t) is a set of some of the vertices that can be reached in MGnp 
from the set .4,(0) by a path of length t or less. It is possible that an inactive 
set Si{t) becomes empty, which means that in MGnp there are no loops at 
the corresponding vertex i, and all d edges incident to i are closed. 

Given the information on the active points and their partners chosen and 
deleted in the first t steps, the random pairing on the set of nd — 2t remaining 
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points remains uniformly distributed. In other words, s" is chosen uniformly 
at random in the set UiG[n] ^ii^) \ {^'l- 
Introduce 

Ij{t) = \{i^l{t):\Si{t)\=j}\, 0<j<d, 

the total number of currently inactive vertices i such that the points from 
Si = Si{0) have been chosen d — j times as partners s" of active points s' 
during the first t steps. [The reason that those z's are still inactive after 
t steps is that the corresponding d — j pairs {s',s") have all been closed.] 
Introduce also 

j 

the total count of currently inactive points, and 

A{t)= Y: \Sr{t)l 

the total number of currently present active points; in particular, ^(0) = 
d\A{0)\. Let = {Ij(t)}o<j<d. It is easy to see that {X(t)} = {A{t),I{t)} 
is a Markov chain. Assuming A{t) > 0, let us compute the one-step transition 
probabilities. There are three kinds of transitions: 

(a) s" is currently active. The (conditional) probability of this transition 

is 

A{t)-1 + I{t)' 

The next state is 

A{t + l)=A{t)-2, I(t + l)=I(t). 

(b) s" is currently inactive, that is, s" G Si{t), i ^Z{t), \Si{t)\ = j, for 
some j G [l,d]. There are two alternatives. 

(bl) (s', s") is open. The probability of this transition is 

jlj{t)p 
A{t) - 1 + I{t) ■ 

The next state is 

A{t + l)=A{t)+j-2, 
h{t + l)=h(t), £^j, 
I,(t + l)=/,(t)-l. 
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(b2) (s', s") is closed. The probability of this transition is 

jlj(.t)q 
A{t)-1 + I{t)' 

The next state is 

A{t + l)=A{t)-l, h{t + l)=h{t), i^j-l,j, 

+ 1) = + 1, Ij{t + 1) = Ij{t) - 1. 

In all cases, 

(3.1) A{t + l) + I{t + l) = A{t) + I{t)-2, 

as it should be, of course. If A{t) = 0, then X(t + 1) = X(t), that is, every 
(0,1) is an absorbing state. 

Note. This Markov chain is a close relative of the Markov chain for the 
bootstrap site percolation in [4]. There are substantial differences though. 
One is that in [4] |^(0)|, the number of active vertices at the start, was of 
order n, while here we will have to consider |^(0)| = 1 or 2, that is, to start 
with just one or two active vertices. With |^(0)| that small, the process runs 
out of active vertices relatively soon, with probability P„ bounded away from 
zero, or even approaching 1, if n^/^(p — p*) -/-^ oo. Since there are n vertices, 
our task basically is to find a way to handle something like n(l — P„) for 
various ranges of p. 

To continue, we average over the three possibilities, (a), (bl) and (b2), 
and obtain the equations for the conditional expectations E[X(t+ l)|X(i)] = 
E[X(t + l)|o]: if A(t) >0, then 

. OA 



E[A(t + l)|o] =A(t) + 



A{t) - 1 + l(t) 



(,+ 1)1,^1 (t)g 

Id+i(t) = 0. + 1) = /j(i) — 1 when s" belongs to an inactive set of car- 
dinality j, no matter whether {s',s") is open or closed.] We will denote the 
right-hand side of (3.2) by X(t) + R(X(t)). 
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3.2. Differential equations approximation. 

3.2(a). As we remarked earlier (note following Lemma 1), whp the num- 
ber of isolated vertices in MGnp is of order n exactly. This means that whp 
I{t) is of order n for every t, suggesting to scale X(t) by n. Another argu- 
ment in favor of such scaling is that, if not for —1 in A{t) — 1, -R(X) would 
have been homogeneous, zero degree, vector-function of X. Besides 

ft + 1 t 
l = {t + l) -t = n{ 

\ n n 

Hoping that, when it matters, X(f -|- 1) is relatively close to its conditional 
expectation E[X(i -|- l)|o], (3.2) leads us to conjecture that the random se- 
quence {X(t)} is well approximated by {nx(t/n)}. Here x(r) = (a(r),i(r)) 
is the solution of the (deterministic) differential equations 



a 



' - ^-2)+p^^(j-2)+g5:-^(-i), 



(3.3) 

,' = ^Al±}h±l^, 0<j<d, 
■' a + i 

where Zd+i(r) =0, subject to initial conditions 

(3.4) a(0) = ^, ij{0) = 0, 0<j<d-l, i^{0) = l-^. 

n na 

We will denote the right-hand side of (3.3) by Ro(x). 
Let us solve (3.3). It follows from the equations that 

(a(r)+z(r))' = -2; 

in view of (3.1) this is hardly surprising. As o(0) +i{0) = d, we have 

a(r) + i{T) =d-2T. 

Therefore we consider r < d/2 only. Introducing u = ln(a(r) -|-i(r))~"^/^, and 
fj{u) := ij{T), we obtain a system of linear (birth-and-death type) equations 

du 

By a backward induction on j it follows that, for all ui, U2, 



= lU + - jfj , 0<j<d, fd+i {u) = 0. 



(3.5) f,{u2) = e-^" 1^"' ■ (1 - e~^Y~^fr{uil u = U2-u^. 

r=j ^ 

Or, setting f (•) = (/o(-), • • • , /d(-)), 

(3.6) f(n2) = M{unui), Mjriu) = e-^W' (^) " e"")'"'l{r>,}. 
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It follows that 

(3.7) M{vi + V2) = M{vi)M{v2), M {-v) = M {v)~^ . 

For lis = ln(a(r) +i(r))-i/2, m =ln{a{0) + iiO))-^/"^ , (3.5) yields 

d 



r=i 

(3.8) 



f a{r)+i{T) \ V2 ^ / 2r \ 

U(0)+^(0); ^ a{0)+mJ ■ 

For the initial conditions (3.4), the formula (3.8) becomes 

ij{T) = djy{Ty[q{l - y{T))f'Ua{0), 0<j<d, 



(3.9) 



Therefore 

(3.10) i{T) =Y.Mr) = Mr){py{r) + qf-^idi^), 

3 

and 



(3.11) 



a(r) = (a(r) + z(r)) - i(T) = S{y{T),p), 
S{y{T),p) := (iy(r)[y(r) - (py(r) + g)"'"^id(0)]. 



We will need to consider ^(0) equal to d or 2d only, corresponding to a 
single active vertex, or to two active vertices at t = 0. So in (3.10)-(3.11), 
^d(O) = 1 — n^^ or id{0) = 1 — 2n~^. S{y,p) = has two roots, and y = 
y{p) G (0, 1). Intuitively we anticipate that the random process will run out 
of active points, thus will terminate, at a moment close to nf, 

(3.12) f=T{p) = ^ii-fy, 

see (3.11), second line, f = 0(l/n) for p <p* , and f is bounded away from 
zero if p> p*. Our task is to find a rigorous argument that will also cover 
the case p — p* ^0. 

3.2(b). For U2 = 0,ui = ln(a(r) + i{T))-'^/'^, (3.5) becomes 
F,(x(r))=i,(0), 

(3.13) F,(x) :=y-^J2<l'-^ (i) - y~'y~''r, 0<j<d, 

r=j ^"'^ 

y = (a + i)i/2_ 
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In a vector-matrix form, the second line in (3.13) can be expressed as 

(3.14) F(x) = M(ln(a + i)^/2)i, x=(a,i). 

While (3.8) solves the initial value problem, (3.13) provides a collection of 
d+1 functions of x, which are the integrals of (3.3). That is, we have 

(3.15) grad*Fj(x) • Ro(x) ee 0. 

Next we use the integrals Fj{x) and (3.15) to construct a collection of expo- 
nential supermartingales, which will allow us to bound probability of "large" 
deviations of X(t) from nx(t/n). 

Lemma 2. Let 7 = j{n) 00 and 7 = 0{n). Introduce 

Qj(t) =exp{7[F,(?i-iX(t)) - F,(n-iX(0))]}, < j < d. 
There exists an absolute constant c > such that 

]E[Qi(t+l)|o]<(l+c(7/n)2)Q,-(t), 
whenever I{t) > ne^ jl. 

Note. I(t) is bounded below by the number of isolated vertices in 
MGnp, which is asymptotic, in probability, to nq'^. So whp the condition 
in the lemma holds for all t. 

Proof of Lemma 2. Obviously we need to consider the case A{t) > 
only. -Fj(x) is certainly twice continuously differentiable for a + i> (f' jl. 
Therefore, since X(t + 1) — X(t) is uniformly bounded, 

= exp{7[F,(n-iX(t + 1)) - F,{n-^Xm\ 

= exp{7[(gradFj(n-^X(t)))*(n~^X(t + l) -n~^X(t))] 

+ 0(7n-2)} 

= 1 + 7gradFj(n-^X(t))*(n-^X(t + 1) - n-^X(t)) + Oiri^ jn^). 
Notice that 

E[(gradF,(n-^X(t)))*(?i-^X(t + 1) - n-^X(t))|o] 

= (gradFj(n-^X(t)))*E[n"iX(t+ 1) -n-iX(t)|o] 
= (gradF,-(n"iX(t)))*n"iR(X(t)) 
= (gradF,(n"iX(t)))*n^iRo(n"iX(t)) +0(n-2) 
= 0{n-\ 
see (3.15). Consequently 

E[g^-(t + l)|o] =gj.(t)[l + 0(7/n2 + 72/^2)] <Q^-(t)(l + c72V). □ 
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Lemma 3. Let T be the first time t when either A{t) =0 or I{t) < nq'^/2. 
(i) Then for any fixed s, and z > 0, 

(3.16) p|max|Fj(n-iX(t)) -Fj(n"^X(0))| > z| < 26"^"+'^^^'/"', 
Ts=T As. (ii) Consequently, for some absolute constant c\ > 0, 

(3.17) ^jmax \n-^^{t) - x(n-4)| > z| < 2{d + i)e-^^-i^+^^^Vn\ 

Note. For (3.17) to be of any help, 72; should substantially exceed 
s7^/n^, and moreover will have to grow logarithmically with n. On the 
other hand z, needs to be much smaller than x{n"^t), for i's in question, 
otherwise (3.17) is pointless. Selecting 7, z, which meet these conflicting 
requirements, will depend on the range of p under consideration. 

Proof of Lemma 3. (i) Define 

Qj{tATs) 



Qj{t) 



(1 + c72n-2)tAT. • 
By Lemma 2, {Qj{t)} is a supermartingale, that is, 

E[gj(t + l)|o] <Qj-(t), t>0. 
Let Tz be the first time t<Ts such that 

Fj{n"^y^{t)) - Fj(n^iX(0)) > z, 

and set 7^ = + 1 if no such t exists. Now 7^ is a stopping time. So, using 
the Optional Sampling Theorem (Durrett [14], Chapter 4), 

E[Q,(T,)]<Q,(0) = 1. 
On the event {Tz < T^} we have 



So 



(1 + C7"'n "')** 



p|max(F,(n-iX(t)) - F,(n-iX(0))) > z| 

<P{Q,(T,) >e^^-^^^''^"'} 
E[Q,(T,)] .,_,,^2„-2 
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In exactly the same way we obtain 



p|max(-F,(n-^X(t)) + F,{n~'X{0))) >z\<e 



So (3.16) follows. 

(ii) According to the definition of Fj^s, we have thus proved that 



(3.18) 



p|max||M(C/(t))(n-il(t))-Af(C/(0))(n-il(0))||<z 
>l-((i+l)e" 



(||f|| ■.= maxj\fj\). Here U{t) =ln{n-^{A{t) + I{t))y/'^ and M(-) is given by 
(3.6). On the event £"„ in (3.18) we have: for t <Ts, 

M{U{t)){n~^I{t)) - M(C/(0))(n~^I(0)) = Z, ||Z|| < z. 

Now M{U{t)y'^ =M{-U{t)), and 

U{t) < ln{n-\A{0) + /(O)))^/^ = Ind^/^ ^ 

Using the definition of M(-) in (3.6), we see then that ||M(— ?7(t))|| < m, for 
some absolute constant m > 0. So on E^, 

mz > ||n"^I(t) - Ar\U{t))M{U{0))n~^I{0)\\ 
= ||n~^I(t) - M{U{0) - U{t))n~^I{0)\\. 

Here 

rr^r.^ rr,^ , , n"M(0) + n-l/(0) \ 2t \ "^^^ 



so 

M(C/(0) - [/(t))n-^I(0) = i(t/n); 
see (3.8), (3.9). Therefore it follows from (3.18) that 

pj^rnax |[?i^^I(t) — i(t/n)|| < mz^ 

This together with 

n~'^A{t) + n~^I{t) = a{t/n) + i{t/n) = d- 2t/n 
imply (3.17). □ 



EDGE PERCOLATION ON A RANDOM REGULAR GRAPH OF LOW DEGREE 23 

4. Proof of Theorem 1 (Subcritical case). Let C„ be size of the com- 
ponent of MGnp which contains a given vertex, 1 say. This component can 
be determined through the percolation process with the starting active set 
^(0) = {1}. In Lemma 3 we introduced the stopping time T, the first t such 
that either A{t) (the number of currently active points) is zero, or I{t) (the 
number of currently inactive points) is below nq'^/2. Clearly C„ < 1 + T. Set 
in Lemma 3 

'y = a{p*—p)n, s = (3{p* — p)~^lnn, 
a > 0, /3 > to be specified shortly. Pick z > such that 

(4.1) 7z>2cs^ z>2ca— — 

in which case 

7' 2 
jz — cs—^ > ca pinn. 

Then, by Lemma 3, 

(4.2) pjmax |n-^X(i) - ^(n-^t)\ > z j = ©(n"^""^) = oijT^), 

if > 1, the restriction we meet by selecting /3 > c^^a"^. On the event 
En\^ complementary to that in (4.2), we have 

(4.3) \n-^A{J,)-a{T,ln)\<z. 
Here, by (3.11), 

a(r) = (iy(r)[y(r) - {pyir) + qf'^il - n~^)], 

(4.4) 

y(r) = (l-2r/d)V2. 

Suppose that, on -Eni) T > s. Then (4.3) implies that z + a{s/n) must be 
positive. Since 

(4 5) /3 ^^'^ (3 ^'^'^ 

n n{p*—pY ni/3(ni/3(p* — p))2 ' 

as v}/'^(p* — p) — > 00, we have 

1 - y{s/n) = 1 - (1 - 2s/{nd)f''^ ~ s/{nd). 

Therefore, by convexity of y — {py + q)'^~^{l — n"^), 

a{s/n) < dy{s/n){n~^ + [I - p{d - - n~^)]{y{s/n) - 1)} 

= -dy{s/n)[{p* - p)(l - y{s/n)) + 0{n~^)] 

S{P*-P) 



n 
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as (p* — p)s ^ oo. So z + a{s/n) < provided that 
fA ft^ ^ ip*-p)s 

(4.6) 

We meet the restrictions (4.1) and (4.6) by choosing, for instance, 

^ = 0.5(2ca + 0.5)s(p* -p). 

With z so defined we come to the conclusion that T < s on the event Eni ■ 
Thus 

F{Cn >l + s)<P{T>s)< P(^^i) = o(n-i). 

Using the union bound we see that whp MGnp does not have a component 
of size exceeding 1 + P{p* — p)"^ logn, that is, 

(4.7) < l+/?(p* -p)-2lnn. 

This completes the proof of Theorem 1. 

Note. Observe that s/n^O if n^/'^{p* — p) » Inn; see (4.5). That is, 
the statement holds under this weaker condition. However, Theorem 2 makes 
it clear that when p* — p is of order n~", a G (1/3, 1/2), the upper bound for 
Ln^ is roughly n?/^ , thus smaller than our current s by n^^'^"^/^) factor. 

5. Proof of Theorem 2 (Nearcritical case). Upper bound. The argument 
parallels the previous proof. This time we set 

(5.1) 7 = n2/3(lnn)", s = n^/'\\nnf, 

for a > 0, /? > to be specified. The counterpart of (4.1) is the condition 

, , S7 (lnn)"+'3 
5.2 z>2c^ ^ z>2c^ 

From (4.4) it follows that 

a{s/n) =dy{s/n) n"^ + (1 -p(d- 1)(1 - n~^))(y(s/n) - 1) 

-P ("^ 2 ^) (1 - n'^){y{s/n) - if + 0((1 - y{s/n)f) 

^-dpd- I2{y{s/n) - if ~ -d-^p (^^ 2 ^ ) 

provided that s/n ^ |p* — p|. Since |p* — p| < n~^/'^{hin)'^ , we satisfy this 
condition by choosing in (5.1) (3 > a. Then z + a{s/n) < if 

(5.3) ,<i^d-l){d-2)'- = {d-l){d- 2f-^. 



'd-l\ ^ 
n^ 
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Comparing (5.2) and (5.3) we see that 

(lnn)"+^ 
^"^"^ n2/3 

satisfies both conditions if a < /3. For this choice of z, 

if 2a + /3 > 1. If /? > max{l/3, a}, we make 2a + /? > 1 by selecting a suffi- 
ciently close to P from below. For these a, /3, we have thus proved that 

lim P(lW > 1 + n^/^ilnnf) = 0. 

6. Proof of Theorem 3 (Supercritical case). Now 

(6.1) uj = uj{n) ■.= n^/^{p-p*) > (Inn)'^^"), 

where a(n) — > oo however slowly. 

Preliminaries. We will have to consider two cases, with one and two 
initially active vertices. So this time 

a(T) =dy{T)[y{T) - {py{T) + qf~^{l - m/n)], 

(6.2) 

y(T) = (l-2r/d)i/2, 

where m = 1, or m = 2. 

Let y = y{p) be the root of fm/niv) = 0, 

fm/niy) ■=y- {py + qf'^i'^ - m/n), 

so that 

is the root of a(r) = in (0,d/2). Since fm/n{y) ^ /o(y) and 

(6.3) -{d-m- 2)/ < < -\[d -l){d- 2)p\\ 
it follows easily that 

(6.4) \y-Ti\=0(n-^l'^), 

where vr G (0, 1) is the root of /o(y) = in (0, 1). From (6.3) for m = 0, it 
follows directly that 1 — vr is of order p — p* exactly, which we express by 
writing 
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Then, by (6.4), 

(6.5) fxp-p*, 

too, as p — p* ^ n~^^'^. We will also need bounds for a(r) when r is (rela- 
tively) close to or f . The convex function fm/niu) attains its maximum at 
a point y € (y, 1), and a simple calculus shows that 

l-y 1 

- = - + 0{p-p*) =^ y-y>:p-p*. 

l-y 2 

Using this fact, /^/„(y) = 0, and (6.3), we obtain that 

\y — y 

I [P-P ), 

f'm/niy) 




p — 

l-y ^ 

p — p* 

if eo > is sufficiently small. Therefore, as fm/niv) = fm/n{l) = 0; 

fm/n{y)>^-{y -y){p-p*), — — T<eo, 

p — p* 

fm./n{y)'^{i-y){p-p*). — T-'^o- 

' p — p* 

Since a{t) = dy{T)f.yn/n{y{'T))^ a little reflection shows then that 
min{a(r) : r € [e(p — p*), t — e(p — p*)]} > ce(p — p*)^, 

(6.6) 

a(f + e{p-p*)) < -ce{p - p*)'^ , 
if e > 0, c > are sufficiently small. 

Running time dichotomy. Let m = l. In Lemma 3 set 

(6.7) 7 = n2/3, s = rif(l + ;3Lj-^), 

/? > being fixed. The analogue of the conditions (4.1), (5.2) is 

7z>2cs^ z>2cf(l + (3uj~^)n~^/^, 

which is met (for n large enough) if 

(6.8) z > cujn~'^l^ 

[see (6.5)], where c' is independent of /?. Further, using (6.5) and the second 
line in (6.6), 

a{sln) = a(f + Pfuj-^) < -c"^. 
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So, picking /5 > 1 sufficiently large and setting z = fi^^'^ujn"'^/^ , say, we ensure 
that z satisfies (6.8) and 

(6.9) z + a{s/n)<Q. 

Since z = z{j3) satisfies (6.8), by Lemma 3(ii) we obtain: there exist Cj = 
Cj(/?) > 0, such that 

p{||n-ix(r,)-x(n-ir,)|| >z} 

(6.10) > 1 - exp(-ci7z) 

> 1 - exp(-C2Cj) = 1 - 0{n^^) > 0. 

Let us denote the event in (6.10) by On S„ n {T > s} 

z + a{s/n)>n~^A{Ts)>0, 
which contradicts (6.9). Hence 

Bn G {T < s} n {|n~^X(r) - x(n~^r)| < z}. 
In particular, since A(T) = 0, 

(6.11) z > \n~'^A(T) - a{n-^T)\ = \a{n-^T)\. 
Suppose that 

— (P-P*) < — <f-—{p-p*). 
coo n CLO 

Then, by the first line in (6.6), 



oj n 



2/3' 



and, by (6.11) and z = P^^'^ton we must have > (3, which is impos- 
sible since (3 > 1. Therefore 

Bn <^{T<s}n {n'^T ^ [c-^Pn-^^^,f - c'^ f3n-^/^]}, 

as p — p* = n~^^^uj. Or, recalling that s = nf{l + /3u;~^), f ^p — p* , 

(6.12) Bn C {n^^T € [f + an"^/^,f - an"^/^]} U {n~^T < an^^^^j, 

for some constant a > 0. In words, with probability superpolynomially close 
to 1, the random percolation — triggered by a fixed initially active vertex — 
either stops after 0(n^/^) steps, or runs much longer, for about nf x wn^/^ 
steps. 
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... and uniqueness of a giant component. On 

:= Bn n {n~^r G [f + an'^/^, f - an'^/^]} 

we have [see (3.9)] 

j=0 j=0 

= (pyin-^T) + - n~^) + 0(u;?i-2/3) 

Therefore C„, the size of the component that contains our initially active 
vertex, on B* satisfies 

d 

Cn = n-Y, Ij{T)n = n[l - (pvr + g)'^] + 0{n'''') 

j=0 

(6.13) 

X Lor?!"^ = n{p — p*). 

And, of course, C„<1 + T<1 + an^/^ on 5„ n {T < an^/S}. Thus, with 
probability 1 — 0{n~^), \/K > 0, either |T — nf| < an^/^ and C„ ~ — 
{p'n + qY],OTT<l + av?/'^ . 

We observe now that no changes whatsoever are needed to show validity 
of this claim for T and C„, the running time and the terminal number of 
active vertices of the percolation process started by two fixed vertices u and 
V. That is, with probability 1 — 0{n~^), V-ftT > 0, 

(6.14) C'„<l.ln[l-(p7r + (/)'^]. 

And the event (6.14) excludes a possibility that u and v belong to two 
different components, each of size relatively close to n[l — (pir + q)'^]. Picking 
K > 2, we obtain then: with probability 1 — 0{n"^^~^^), there are no two 
vertices u,v [n] such that the two components containing them are disjoint, 
and both are of size of order n{p — p*). Putting it differently, whp there may 
exist at most one "giant" component, and all other components have size 
0{n'^/^) at most. 

Other components. . . . Let us show whp that nongiant components are 
actually much smaller, of order n'^^^uj~'^{\nny^°^^\ 
To this end, define a sequence {7^, Sfc, Zfc}. For k>0, 

and 

(6.15) Zk > 2cso^ = 2caa;'^'=-^^ 
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fjfc > 0, /3fc > 0. In particular, Pq = 0, so that sq = ar?!'^ is the (probabihs- 
tic) upper bound for ah nongiant components already proved. Further, 
increases with fe, and 

(6.16) c7fc = — + — j , p>\. 

I ma; 

So 

IkZk - csfe4 = cacu^-^-^* = (Inn)''. 
Applying Lemma 3 with 7^, and Zfc, we obtain that 

(6.17) |n-iX(t)-x(n"it)| <Zfc Vt < r,,(= min{r, Sfe}), 
with probability 1 - 0(e-(i°")''). Now 

^ = 0^-1^0; 

p — p* 

so, using (6.6), 

(6.18) min{a(t) :t G [sfc+i/n, Sfc/n]} > a(n~^Sfc+i) 
Pick /)' S (l,p), and set 

^« 1 o^ R 1 _L In Inn 

(6-19) /3fc+i = l + - . 

z zma; 

Then [see (6.16)], 

/5fc+i - 1 - (/3fc - o-fe) = -^^^ — , 

imu; 

which means that we can pick 2;^ satisfying (6.15), and such that 

Zfc = o(a(n"^Sfc+i)); 

see (6.18). But then [sfc+i,Sfc] on the event (6.17), since otherwise, by 
(6.17), we obtain that 

a(n-isfc+i) > |a(n-i(T))| = \a{n-^T) - tT^ A{T)\ < Zk, 

contradicting our choice of z^- 

The solution of (6.19), subject to /?o = 0, is 

p' In Inn . , , / p' In Inn 
Pk = 2- 2-^=+^ 1 ^ 



In a; V 2 In 

For ko = [logglnw], 

(/)' + o(l))lnlnn 



/3feo=2 



Inu) 
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so that 

2/^ n2/3(lnn)^'+°(i) 
Skr,=auj ^''on^''^ = a ^ 4 , 

for all 1 < p' < p. By the union bound inequality, we have then 

P{Te[sko,so]}< ^ P{re 

k=0 

= 0(e-(^'^")''A;o) = 0(e-(i°")'' Inlnn) 
= 0(n-^) VK>0. 

Now So = an^/3 , and we had proved that 

P{T € [an2/3,nf - an2/3]} = O(n--f^). 

Consequently 

P{T € [u;-2n2/3(lnn)i+"(i),nf - an^/^]} = ©(n"^), 

for all K > 0. Here T = r(ii) is the running time for a fixed initial vertex u. 
That is, switching to Cn{u) and using the union bound, 



I we [nl 



(n) - n[l - (pvr + g)^] or = 0(w-2n2/3(lnn)i+°(^))} 



> l-0(n--^+^). 



... and existence of a giant component. Corollary 1 implies that whp 
Gnp* contains a component of size w~^n^/^ at least. As the multigraphs 
Gnp* and Gnp, p> p*, can easily be coupled in such a way that Gnp* C Gnp, 
we conclude that whp Gnp also has such a component. Since 

o;-in2/3» ^-3/2^2/3^ 

it follows from (6.20) that whp this component has to be giant, of size close 
to n[l - [pir + qy]. 
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